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C^ ' Abstract 

o 

oo , 

O , We construct the Lemaitre-Tolman-Bondi (LTB) dust universe whose distance- 

redshift relation is equivalent to that in the concordance A cold dark matter (/ICDM) 

/\ • cosmological model. In our model, the density distribution and velocity field are not 

c^ \ homogeneous, whereas the big-bang time is uniform, which implies that the universe 

is homogeneous at its beginning. We also study the effects of local dumpiness in the 

density distribution as well as the effects of large-scale inhomogeneities on the distance- 

redshift relation, and show that these effects may reduce the amplitude of large-scale 

inhomogeneities necessary for having a distance-redshift relation that is the same as 

that of the concordance /ICDM universe. We also study the temporal variation of 

the cosmological redshift and show that, by the observation of this quantity, we can 

distinguish our LTB universe model from the concordance ylCDM model, even if their 

redshift- distance relations are equivalent to each other. 
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§1. Introduction 

A cold dark matter (ylCDM) models have achieved wide acceptance as concordance mod- 
els due to the results of observations over the last decade. In particular, cosmic microwave 
background (CMB)^-* and supernovae (SNe)^-*"^-* observations have played critical roles in this 
acceptance. The isotropy of our universe is strongly supported by the CMB observations. 
Thus, if we assume that our universe is homogeneous, results from SNe observation suggest 
that the volume expansion of our universe is accelerating. The accelerating expansion of a 
homogeneous and isotropic universe means the existence of exotic energy components, the 
so-called dark energy, within the framework of general relativity (GR). If we consider dark 
energy as a perfect fluid, it has negative pressure. One of the candidates for dark energy is 
the cosmological constant. However, there are several crucial problems regarding the exis- 
tence of the cosmological constant or other dark energy candidates (see, for example, Ref.6)). 
No one knows the origin of dark energy and there has not yet been a conclusive illustration 
of dark energy. 

There are other possibilities that explain the CMB and SNe observations, observations. 
In the above logic, we have imposed two assumptions. One is validity of GR at all cosmo- 
logical scales. The other is the homogeneity of our universe. Hence, in order to describe the 
observational results without dark energy, we need to discard GR or homogeneous models 
of the universe. In this paper, we attempt to describe the observational results without dark 
energy or any modification of gravitational theory, but with inhomogeneous universe models. 

The basic idea is that we are in a large underdense region, i.e., a large void; we reject 
the Copernican principle, which states that we live at a typical position in the universe. 
Pioneering works include those Zehavi et al.^^ in 1998 and Tomita^^"^''-' in 2000 and 2001. 
Zehavi et al. analyze early SNe data and suggested that such a large void might exist 
around us without dark energy. Tomita proposed the void universe model and discussed 
the possibility of explaining the observed magnitude-redshift relations of SNe.^^"^°^ There 
are several works in the same direction^-'^^"-'^^^ (See also the reviews. "'^^-'"^^-'). In these works, 
Lemaitre-Tolman-Bondi (LTB) solutions^^-*"^^-* are often employed. LTB solutions are exact 
solutions of the Einstein equations, which describe the dynamics of a spherically symmetric 
dust fiuid and are useful for constructing universe models with a spherically symmetric void. 

Recently, several authors have discussed the possibility of explaining the CMB observa- 
tions. ^^-'"^'^-' They used asymptotically homogeneous LTB models and reported that these 
models may be consistent with the observed CMB anisotropy. In order to obtain more pre- 
cise predictions on whether we are located near the center of a large void, we need to study 
the evolution of nonspherical density perturbations in the LTB universes. ^^'' Other methods 



of observationally investigating the inhomogeneity of our universe have been proposed in 
several papers. ^^^"^^^ 

The aforementioned works are very important for strengthening the observational foun- 
dation of physical cosmology. At present, the Copernican principle is not based on sufficient 
observational facts. However, by virtue of the recently improved observational technologies, 
we might reach the stage at which we are able to investigate through observation whether our 
location in the universe is unusual or not . In order to test inhomogeneous universe models 
by observations, it is important to know the types of inhomogeneous universe models that 
can be used to explain current observations, and to reveal what predictions are given by 
these universe models. This is the subject of this paper. 

The inverse problem using LTB universe models is a useful method for investigating the 
possibility of explaining observational results with inhomogeneities in the universe. ^^-'' i4)'34),35) 
LTB dust solutions contain three arbitrary functions of the radial coordinate: the mass func- 
tion M, the big-bang time t^ and the curvature function k, and the inverse problem means 
determining M, te and k so that a given distance-redshift relation is realized on the past 
light-cone of an observer. In order to specify the three arbitrary functions, we need three 
conditions. One of these conditions corresponds to the choice of the radial coordinate and 
thus has no physical meaning. Hence, one more condition in addition to the distance-redshift 
relation is necessary. 

Iguchi et al.^^-* showed that the distance-redshift relation in the standard ylCDM model 
can be reproduced in LTB dust universe models. When they solved the inverse problem, 
they imposed two kinds of the additional condition: one is a uniform big-bang time te = 
and the other is a vanishing curvature function A; = 0. It should be noted that their models 
failed to reproduce the distance-redshift relation for z > 1.7. The reason why they could not 
continue their computation beyond z ~ 1.7 is the existence of the "critical point" discussed 
by Vanderveld et.al.^^-* They reported the difficulty in solving the inverse problem and 
concluded that it is unlikely that a solution to the inverse problem can be found in all the 
redshift domain. However, Tanimoto and Nambu recently solved the inverse problem in all 
the redshift domain with a nonuniform t^.^^^ In this paper, we show that it is possible to 
obtain a solution to the inverse problem with the condition of the uniform big-bang time t^. 
*) In terms of perturbation theory, the inhomogeneity of the big-bang time corresponds to 



*■* Mustapha et al.'^^' proposed the following theorem '.Subject to the conditions of Appendix B in Ref.34) 
for any given isotropic observation of apparent luminosity l{z) and number count n{z) with any given source 
evolution function L{z) and total density over source number density rh{z), a set of LTB functions can be 
found to make the LTB observational relations fit the observations. It is nontrivial whether the same 
statement holds for any given set of some quantities different from l{z)^ nlz), L{z) and rh{z). In this paper, 
the big-bang time tB{r) is given as one of the quantities. Thus, the situation is different from the case in 



the decaying mode. The condition of uniform big-bang time might guarantee the consistency 
of the present model with the inflationary scenario, since a universe that experiences inflation 
is almost homogeneous immediately after the inflationary period is over. 

In this paper, we also discuss the effects of local dumpiness. In addition to the effect of the 
large-scale void structure, small-scale dumpiness may affect the observed distance-redshift 
relation. If there are clumpy objects such as galaxies in the foreground of light sources, data 
for the apparent luminosities of these sources are often contaminated with gravitational 
lensing or absorption. For this reason, supernova teams specifically search for objects with 
a minimum amount of intervening material in the foreground. ^^-''^^-' The observed sources 
might be biased by this selection, since the light of the observed sources might propagate 
in a lower-density region of our universe. The simplest way to take account of this effect 
is to introduce the so-called smoothness parameter a, which is the ratio of the smoothly 
distributed matter density except for clumps to the mean energy density p for all matter, 
assuming that the energy density on the paths of observed light rays is given by ap.^^^^^^^ 
Tomita has studied the effects of small-scale dumpiness in a local void model by introducing 
the smoothness parameter a. In this paper, we also investigate the corresponding effects by 
introducing a. 

The time derivative of the cosmological redshift is a very important observable quantity 
in distinguishing LTB models from concordance ylCDM universe models. We study it and 
show that it gives a criterion for observationally deciding whether the universe is described 
by an LTB model with uniform big-bang time. 

This paper is organized as follows. In §2, we briefly review LTB dust universes and list 
basic equations that we have to solve. The temporal variation of the cosmological redshift 
of comoving sources is discussed in §3. Then, we show details of our numerical method, 
particularly, focusing on singular points of the equations at the center and the critical point 
in §4. Numerical results are given in §5. Section6 is devoted to a summary and discussion. 

Throughout this paper, we use the unit of c = G = 1, where c and G are the speed of 
light and the gravitational constant, respectively. 



Ref.34). 



]2. Basic equations 



2.1. LTB dust universe 



LTB solutions are exact solutions to the Einstein equations, which describe the dynamics 
of a spherically symmetric dust fluid and whose line element is written in the form 

ds' = -de + iM^iL^rf^2 ^ R2^t^^)dQ\ (2.1) 

1 — k[r)r^ 

where k{r) is an arbitrary function of the radial coordinate r. The LTB solutions include 
homogeneous and isotropic universes as special cases; in this case, A; is a constant called the 
curvature parameter in the appropriate gauge, and thus we call it the curvature function. 
The stress-energy tensor of the dust is given by 

T'^'^ = p(t, r)M'^M^ (2-2) 

where p(t, r) is the rest mass density of the dust and u^ = Sq is the 4-velocity of a dust 
particle. The Einstein equations lead to the equations for the areal radius R{t, r) and the 
rest mass density p(t,r), 

{dtRf = ~k{ry + ^^ (2-3) 

rt 

and 

where M(r) is an arbitrary function of the radial coordinate r. We assume that p is non- 
negative and that R is monotonic with respect to r, i.e., drR > 0. 

Following Tanimoto and Nambu,^^-* the solution to Eqs. (2-3) and (2-4), which represents 
the expanding universe, is written in the form 

R{t, r) = (6M(r))i/3(t - tB{r)f^''S{x), (2-5) 

where tB{r) is an arbitrary function of the radial coordinate r. The function S{x) is deflned 

as 

cosh ^^—1] — 1 — (sinh ^^—1] — ^—r])^!'^ 

^^^ " 6i/3(sinhv/=^-v/=^)2/3 ' ^ = 6^75 *°^ ^ < °- 

1 — cosw^ [^Jri — 'svvl Jrif'l^ 
'^("^ = 6V3(^-siny^)V3 ' - = W^ f- ->0- (2.7) 
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S{x) is analytic in the domain x < Xc= (tt/S)^'^. Some characteristics of the function S{x) 
are given in Appendix A and Ref.41). We can easily see from the above equations that the 
areal radius R vanishes at t = tB{r). Thus, the function tB{r) is called the big-bang time. 
Following Refs.15) and 20), we define the local Hubble function by 

H{t,r) = ^ (2-8) 

and the density-parameter function of dust as 

In LTB universes, we can define another expansion rate of the spatial length scale, the 
so-called longitudinal expansion rate, by 

H'it,r) = ^. (2.10) 

In the case of homogeneous and isotropic universes, H^ agrees with H. Thus, if we can 
measure the difference between H and H^, it can be used as an indicator of the inhomogeneity 
in the universe. ^^-^ 

2.2. Conditions and equations to determine arbitrary functions 

As shown in the preceding subsection, the LTB solutions have three arbitrary functions: 
k{r), M{r) and tB{r). We have one degree of freedom to rescale the radial coordinate r. To 
fix one of the three functional degrees of freedom corresponds to fixing the gauge freedom of 
this rescaling, and it will be shown later how to fix it. The remaining two functional degrees 
of freedom are fixed by imposing the following physical conditions. 

• Uniform big-bang time te = 0. 

• The angular diameter distance D{z) is equivalent to that in the ylCDM universe in 
all the redshift domain except in the vicinity of the symmetry center in which D is 
appropriately set so that the regularity of the spacetime geometry is guaranteed. 

Here we stress that it is our primary purpose to find an LTB model that fits the observed 
distance-redshift relation well. Thus, it does not matter that the distance-redshift relation 
does not agree with that of ylCDM model only in the vicinity of the symmetry center. 

In order to determine k{r) and M{r) from the above conditions, we consider a past- 
directed outgoing radial null geodesic that emanates from the observer at the center. This 
null geodesic is expressed in the form 

t = t{X), (2-11) 

r = r(A), (2-12) 
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where A is an affine parameter. 

We assume that the observer is always located at the symmetry center r = and observes 
the light ray at t = to- In order to fix the gauge freedom to rescale the radial coordinate r, 
we adopt the light-cone gauge condition that the relation 

t = to-r (2-13) 

is satisfied along the observed light ray. 

Then the basic equations to determine k and M are given as follows: 

1. Null condition 

By virtue of the light-cone gauge condition, the null condition on the observed light 
ray takes the very simple form of 



dr-R = Vl - kr^. (2-14) 



2. Definition of redshift 

The redshift is defined by 



l + z= ^ ■^y^""'-^'^ oc po = -p" = -i, (2-15) 

"M I observer 



where p'^ is the tangent vector of the null geodesic that corresponds to the observed 
light ray and the dot represents differentiation with respect to the affine parameter A. 
By using the freedom to multiply the affine parameter by a constant, we can write, 
without loss of generality, 

t = -r = -i±^, (2-16) 



where 



Ho:=H{to,0), (2-17) 



and we have used the gauge condition (2-13) in the first equality. In this normalization, 
the affine parameter is dimensionless. 

3. Geodesic equation 

One of the geodesic equations for the radial null geodesic is given by 

{drR)i + {dtdrR)i'^ = 0, (2-18) 

where we have used the null condition (2-14). 

4. Dyer-Roeder equation for the angular diameter distance 



As mentioned, in order to fix the remaining functional freedom, we assume the angular 
diameter distance D{z). Then the Dyer-Roeder equation 



dD ^(PD {l + zf 

'^-r + "^ -nr = -47r ^ ^/ apD 2-19 

dz dz^ Hq 



gives us one of the equations to determine the arbitrary functions of LTB universe 
models, where a is the smoothness parameter mentioned in §1. The derivation for 
this equation is given in Appendix B. a is the mass ratio of the smoothly distributed 
components of matter to all the components, and it may vary with time due to the 
formation of structures in the real universe. Therefore, in this paper, we assume that 
a is an input function of the cosmological redshift z. 

Equations (2-14), (2-16), (2-18) and (2-19) are rewritten in the form of five coupled 
ordinary differential equations: 

m = F„(m,A;,r,z,C), (2-20) 

A; = Ffc(m,A;,r,z,C), (2-21) 

r = ^, (2-22) 



c 



dD /dz'' 

-47r(l + zfapD 



HI 



(2-23) 
(2-24) 



where m and C, are defined by 



and 



(2-27) 



6M(r) 
m{r) := ^ (2-25) 

dD 

C ■■= z—, (2-26) 

dz 

respectively. From Eq. (2-4), we have 

_ r^ [3(1 + z)m + HorFrnim, k, r, z, ()] 
^~ 247r(l + z)i?Vl - A;r2 

The derivation of these equations is shown in Appendix C 

§3. Temporal variation of the cosmological redshift 

The temporal variation of the cosmological redshift will give us crucial information about 
the acceleration of cosmic volume expansion or inhomogeneities in our universe. ^^-^ 
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3.1. Homogeneous and isotropic universe 

The line element of homogeneous and isotropic universes is given by 

ds^ = -dt^ + a\t) {dx^ + S{x)dQ^) , (3-1) 

where S{x) = si^X fo^^ ^ closed universe, S{x) = X fo^^ ^ A^-t universe and S{x) = sinhx for 
a open universe. In this subsection, we assume that the universe is filled with dust and dark 
energy characterized by the linear equation of state p = wp, where p is the pressure, p is the 
energy density and w is a constant less than —1/3. We assume that the energy densities of 
both dust and dark energy are nonnegative. 
The Einstein equations lead to 

where oq is the scale factor at t = to, and ^mo and Oxo are the density parameters of the 
dust and the dark energy, respectively. Note that, by assumption, both ^mo aiid Oxo ^^^ 
nonnegative. As is well known, the cosmological redshift z of a light signal emitted from a 
comoving source is given by 

-^-1, (3.3) 

where tc is the time when the light is emitted from the source. The temporal variation of z 
of a comoving source is then given by 

ajto + Ato) a{to) rj f . , , ^A 

Az = — - r^ Ho[l + z- — \ Ato, (3-4) 

a{tc + Atc) a{tc) V Ho J 

where Ate = Ato/ {I + z) and He is the Hubble parameter when the light is emitted from 
the source. Substituting Eq. (3-2) into the above equation, we have 

dz 



,^ Hoil + z) l-^l + nMoz + nxoiil + zy+'^-^-l} . (3-5) 

uto L 

Note that, O^oz is nonnegative, whereas ^xo{(l + 2;)^+^™ — 1} is nonpositive due to the 
assumption w < —1/3. Thus, if the dust is a dominant component of the universe, i.e., 
f^Mo ^ ^XO) then dz/dto is negative. By contrast, in the case of a universe dominated by 
dark energy, i.e., f^xo ^ ^mo; then dz/dto is positive. Therefore, the measurement of the 
temporal variation of the cosmological redshift z will give us crucial knowledge about the 
equation of state if the universe is homogeneous and isotropic. 
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3.2. LTB universe 

In order to obtain the time derivative of the cosmological redshift z in LTB universe 
models, we consider another past-directed outgoing radial null geodesic that is infinitesimally 
close to the null geodesic considered in the preceding section, 



t = tb(A) + 5t(A) and r = rb(A) + 5r(A), 



(3-6) 



where tb(A) and rb(A) denote the null geodesic with the initial condition r(0) = at t(0) = to, 
which was considered in the preceding section. We set the affine parameter so that the 
cosmological redshift is given by 

1 + Z 1 + Zh + Sz 



t = tb + 5t 



(3-7) 



Hq Hq 

where the subscript "b" denotes the value evaluated on the null geodesic (t,r) = (tb,''"b)- 
Thus, we have 

6z = -Ho6i. (3-8) 

Substituting Eq. (3-6) into Eq. (2-18), and taking the first order of 6t and 6r, we have 



6i + 2H^iji + 



Of- OrR . jjL\2 



il6t + 



dtdlR ^^^dlR 



drR 



H' 



Another equation is given by the null condition. 



5f = -5t + H^t^,5t + 



drkr'^ + 2kr dlR 



+ 



drR 



drR 



%5r. 



il5r = 0. (3-9) 



(3-10) 



2{drRy 
The initial conditions for these equations are given by 

(5t(0) = 0, (3-11) 

(5r(0) = 0, (3-12) 

5t(0) = const.. (3-13) 

The temporal variation of the cosmological redshift of a comoving source is given by 

Az{X)=z{X + AX)-zy,{X), (3-14) 

where the infinitesimal quantity AX satisfies the equality 

r(A + Z\A)=rb(A). (3-15) 

From Eqs. (3-6), (3-7), (3-14) and (3-15), we have, up to the first order of AX, 



Az{X) ~ ^^(A) + ib(A)Z\A ~ 6z{X) 



4(A) 
hiX) 



6r{X). 



(3-16) 
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Thus, we have the derivative of the cosmological redshift z with respect to the time of the 
observer as a function of the afiine parameter, 

Since we also obtain z\^[X) simultaneously, we have the relation between z and dz/dtQ. In 
order to solve the differential equations (3-9) and (3-10), we need to know {df.R)h, {d^drR)h 
and {dtd^R)h, which are given in Appendix D. Using Eq.(D-3), Eq.(3-10) can be solved as 

Sr + St = ^ (3-18) 

1 + z 

where we have used Eqs.(2-10) and (2-18). 

§4. Singularity at the center and critical point 

4.1. Resolving singularity at the center 

Suppose that the reference angular diameter distance D{z) is exactly the same as that 
of the concordance ylCDM model with {iluo, ^Ao) = (0.3,0.7), where Oao is the density 
parameter of the cosmological constant. Then, as shown in Appendix E, once i7M(to,0) is 
given, the set of solutions for the differential equations (2-20)-(2-24) is uniquely determined. 
Vanderveld et al.^^-* reported that many of the inhomogeneous models that mimic observa- 
tions of an accelerating universe contain a weak singularity at the symmetry center. This 
singularity is too weak to make the spacetime geodesically incomplete, and thus we may 
accept these models as being effective. 

On the other hand, the accuracy of observations in the low-redshift domain is not suf- 
ficient to uniquely determine the distance-redshift relation. Thus, we need to assume the 
redshift dependence of the angular diameter distance in the low-redshift domain. If the C°° 
model is preferred, we may assume the following input angular diameter distance D{z), 

D = 1^(0.3,0.7) (^) [1 -exp (-2;V52)] +/^(i^^o,o)(^)exp {-z^/5^) , (4-1) 

where 5 is a positive constant, whereas -D(o.3,o.7) and -D(i7mo,o) a-^e the angular diameter dis- 
tances in the isotropic and homogeneous universe with (i^MO; ^Ao) = (0.3, 0.7) and (J7mo? ^ao) 
(^mO) 0), respectively. D{z) is almost the same as the angular diameter distance of the ylCDM 
model for z > 5^ whereas, in the vicinity of the symmetry center, D{z) is almost the same 
as that of the homogeneous and isotropic universe with (^mo, ^aq) = {^mo, 0). It should be 
noted that f]„^o is equal to ^m(^O) 0) and can be regarded as a parameter to specify the LTB 
universe model. 
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Also note that if our vicinity is described well by the very smooth LTB model, the 
distance-redshift relation does not agree with that of the concordance ylCDM model in the 
low-redshift domain. Nambu and Tanimoto have shown that the Maclaurin series oi D{z) for 
the LTB model with a regular symmetry center agrees with the homogeneous and isotropic 
dust-filled universe up to z"^. The angular diameter distance (4.1) is one of the simplest 
assumptions that ensures the regularity at the symmetry center and agrees with that of the 
concordance ylCDM model in the high-redshift domain. 

4.2. Critical point 

The differential equation (2-23) has a singular point at which dD/dz vanishes. Following 
Ref.35), we call this point the critical point and denote the cosmological redshift at the 
critical point by z^r- We can require that the value of ( should vanish at the critical point so 
that the solution is regular at this point. This gives a constraint on the free parameter ^mo. 

The symmetry center r = is another regular singular point of the differential equations 
(2-20)-(2-24). The Runge-Kutta method, which we have used, is generally unstable for 
solving ordinary differential equations toward a regular singular point from a regular point, 
and thus we start the numerical integration from these regular singular points; we numerically 
integrate the equations from the symmetry center r = but not to the critical point z = Zcr, 
and we also integrate them from the critical point z = z^x but not to the symmetry center 
r = 0. Instead, specifying f^mo and the values of m, k and r at the critical point z = z^r, we 
numerically integrate the differential equations outward from the symmetry center r = and 
inward from the critical point z = Zcr to the matching point z = z-^a located in the domain 
between these singular points. 

If we fail to choose appropriate values of i7jno and values of m, k and r at the critical 
point z = Zcr, the resultant solutions are discontinuous at z = Zj^. Thus, we have to search 
for the appropriate initial values for ^mO) and m, k and r at the critical point z = z^x so 
that the following matching conditions are satisfied: 

fc|2 = 2m + = ^|2 = 2m-05 / . „n 

'"|2 = 2m + ^ ''"|2 = 2m-0 5 

S|2 = 2m + ^|2 = 2m — 0- 

Note that if the above conditions hold, the smoothness of the solutions is also guaranteed, 
since the equations for these functions are first-order differential equations. 

We have searched for appropriate initial conditions that guarantee the matching condi- 
tions (4-2) by using the four- dimensional Newton-Raphson method. Using this procedure, 
we can uniquely obtain the solution if the value of 5 in Eq. (4-1) is fixed. 
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§5. Numerical results 

We have solved Eqs. (2-20)-(2-24) by using the numerical procedure described in the 
previous section. In the following subsections, we express k{r{z), m{r{z)), Quito, r{z)), 
H(to,r{z)), H^{tQ,r{z)) and p(to,r{z)) as functions of the cosmological redshift z. We 
also express the time derivative of the cosmological redshift dz/dto as a function of the 
cosmological redshift z itself together with that of the concordance ylCDM universe with 
(f2Mo,^Ao) = (0.3,0.7). 

5.1. Results without local dumpiness (a = 1) 

In this subsection, we assume a = 1. As mentioned above, the solution is uniquely given 
by our numerical procedure if the value of 6 in Eq. (4-1) is fixed. First, we show m{r{z)) and 
k{r{z)) as functions of z in Figs. 1 and 2, respectively. The results do not strongly depend 
on the value of 6 except when z < 6. 
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Fig. 1. ra(r{zy) depicted as functions of the cosmological redshift for various values of b. 
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Fig. 2. k{r{z)) depicted as functions of the cosmological redshift for various values of 6. 

In order to determine the physical properties of the solution, we depict f]M(to,r{z)), 

p(tQ,r{z)), and H{tQ,r{z)) and H^{to,r{z)), respectively, as functions of z in Figs. 3-5. We 

can see from these figures that the resultant inhomogeneity is a large-scale void structure. 

We note that H{to,r{z)) = H^{to,r{z)) in homogeneous and isotropic universes, whereas 
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Fig. 3. Density-parameter functions f2M{to,r{z)) depicted as functions of the cosmological redshift 
for various values of 6. 




Fig. 4. Rest-mass densities p{to, r{z)) depicted as functions of the cosmological redshift for various 
values of 6. 
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Fig. 5. Local Hubble function H{to,r{z)) and longitudinal expansion rate H^{to,r{z)) depicted 
as functions of the cosmological redshift for 6 = 0. 



H{to,r{z)) and H^{tQ,r{z)) are different from each other by about 10% for 2 < 2; < 10 
in the inhomogeneous case depicted in Fig. 5. This result means that, in order to fit the 
distance-redshift relation of the LTB model with observations that almost agree with that 



14 



predicted by the concordance ylCDM universe with (^mo? ^ao) = (0.3, 0.7), the scale of the 
inhomogeneity should be at least a few Gpc. 

5.2. Results with local dumpiness 

In this subsection, we fix the value of 5 as and study the effects of local dumpiness. The 
vanishing 6 implies that D{z) agrees with the angular diameter distance of the concordance 
ylCDM model, and thus the weak singularity appears at the symmetry center. 

In the early universe, the matter distribution might have been smooth, whereas it might 
have been highly clumpy after the structures formed. The simplest way to take account of 
the effects of local dumpiness and the growth of structures is to introduce a z-dependent 
smoothness parameter a{z), which represents the fraction of matter spreading out almost 
homogeneously. The light rays propagating from standard candles to us pass through regions 
filled with matter of energy density a{z)p. 

At present, the 2;-dependence of a{z) is still unclear. However, as the structures grow, the 
spreading-out component of the mass might decrease, and therefore a might be an increasing 
function of z. Furthermore, there might be a typical redshift z = [3 at which almost all of 
the mass components form clumpy structures. This typical redshift f3 will depend on the 
scenario of structure formation. For z < j3, the light rays might propagate to us through 
almost empty regions. Thus, we assume the following form for the smoothness parameter 
a{z): 

a{z) = l-exp[-z^/p'^]. (5-1) 

A bundle of light rays propagating through a region where a = 1 is called a filled beam, 
while a bundle of light rays propagating through a region where a = is called an empty 
beam. Thus, all bundles of light rays are empty beams in the model oi f3 = oo. In the case 
of the empty beam, the angular diameter distance does not have a maximal value, and it is 
clear that the distance in the /3 = oo case cannot fit the distance in the ylCDM universe in 
all the redshift domain. The results depend on the value of f3 as shown in Figs. 6—9. The 
resultant inhomogeneity is also a large-scale void structure. 

We show that it is possible to reduce the amplitude of the large-scale inhomogeneity by 
choosing an appropriate value of /3. In Fig. 10, values of H(to,r{z)) and H^{to,r{z)) are 
depicted as functions of z ioi {3 = 1.1. In this case, the difference between H(to,r{z)) and 
H^{to, r{z)) is a few percent when z '>2, and the amplitude of the inhomogeneity is smaller 
than that depicted in Fig. 5. This means that the appropriate value of the redshift that 
characterizes the period of local dumpiness formation may suppress the size of the void by 
a few Gpc. In our model, this appropriate value of the redshift is given by /9 ~ 1.1. In the 
region 2; > 4 of our model with (3 = 1.1, the geometry of the universe is almost the same as 
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Fig. 6. Rest-mass densities p{tQ,r{z)) depicted as functions of the cosmological redshift for 6 = 
and various values of (3. 
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Fig. 7. Density-parameter functions f^Mitoii^iz)) depicted as functions of the cosmological redshift 
for 5 = and various values of (3. 

that of the isotropic and homogeneous universe. 

5.3. Time variation of the redshift 

We have shown in the preceding subsections that it is possible to construct the LTB 
universe model with the same distance-redshift relation as that of the concordance /ICDM 
model. Thus, it is very important to study how to observationally distinguish these two 
models from each other. Here, we show that the temporal variation of the cosmological 
redshift is a useful observational quantity. In Fig. 11, we depict the derivative dz/dto of 
the cosmological redshift with respect to the time to of the observer at the symmetry center 
r = as a function of the cosmological redshift itself. 

As can be seen from this figure, dz/dto is positive for < 2; < 2 in the case of the 
concordance /ICDM model, while it is negative for all z in the LTB universe models with 
the uniform big-bang time. Therefore, if we observe whether dz/dto is positive or negative 
for 2; < 2, we can distinguish our LTB model from the concordance ylCDM model. From Eq. 
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Fig. 8. Local Hubble function H{tQ,r{z)) depicted as functions of the cosmological redshift for 
5 = and various values of /3. 




Redshift 



Fig. 9. Longitudinal expansion rate H^{to,r{z)) is depicted as functions of the cosmological red- 
shift for (5 = and various values of /3. 



(3-5), we have (i2;/(ito|z=i ~ 0.24:Hq for the concordance ylCDM model, and thus, the variation 
of the cosmological redshift in one year is Az\z=i ~ 1.8x 10^^^(/7o/75 km/s/Mpc). Thus, over 
ten years, Az\z=i is larger than 10~^° for the concordance /ICDM universe model, and this 
value will become observable in the near future as a result of technological innovations.^^-'' '*^'* 

§6. Summary and discussion 

In this paper, we have attempted to solve the inverse problem to construct an LTB 
universe model that has the same distance-redshift relation as that of the concordance /ICDM 
model with (i?MO) ■^^yio) = (0.3,0.7), and we obtained solutions by numerical integration. In 
the present study, assuming an inflationary period in the early universe, we have restricted 
ourselves to models with uniform big-bang time, and the resultant universe model has a very 
large void whose symmetry center is at the observer's position. We have also studied the 
effects of local dumpiness by introducing the smoothness parameter a, and have shown that 
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Fig. 10. Local Hubble function H{tQ,r{z)) and longitudinal expansion rate H {tQ,r{z)) depicted 
as functions of the cosmological redshift for 6 = and /? = 1.1. 
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Fig. 11. Time derivatives of the cosmological redshift for the various models depicted as functions 
of the cosmological redshift. 



the local dumpiness may reduce the amplitude of the large-scale inhomogeneity of the void. 
Our results imply that it is possible to construct an inhomogeneous but isotropic universe 
model with a distance-redshift relation that agrees quite well with the observational data of 
the distance-redshift relation. 

Our LTB universe model is regarded as an unnatural model from the viewpoint of the 
Copernican principle, because the observer stands exactly at the center of the isotropic uni- 
verse. The extent to which we can separate ourselves from the center and remain consistent 
with current observations is discussed in Refs.8) and 45)-48). The strictest limit is given by 
Alnes and Amarzguioui in Ref.45) as 15 Mpc from the observation of CMB. This value is 
much smaller than the cosmological scale, and thus we should remain at a special position if 
the void universe is real. However, no observational data has yet been reported that entirely 
excludes inhomogeneous universe models. Therefore, it is important to know the type of 



inhomogeneous universes that can explain current observations and to propose observational 

methods for testing inhomogeneous universe models. 

In this paper, we have also studied the temporal variation of the distance-redshift relation 

in our LTB universe model whose distance-redshift relation is the same as the concordance 

ylCDM model with [f^uo, ^ao) = (0.3, 0.7). The result implies that if we can observe the time 

derivative of the cosmological redshift with sufficient accuracy, we can distinguish our LTB 

model from the concordance ylCDM universe model. Innovations in observational technology 

might provide us with data on the time derivative of the cosmological redshift in the near 
future.3i)'43),44) 

Finally, we should note that if we ignore the inflationary paradigm, the functional freedom 
of the big-bang time tB(^) is returned, and we might be able to construct an LTB universe 
model with the same redshift-distance relation and, furthermore, the same dz/dt^ as those of 
the concordance ylCDM model by choosing an appropriate big-bang time tB(^).^^'' However, 
in this case, we may need an other mechanism to explain the results in CMB observations 
and other cosmological problems than the standard cosmology starting from the inflation. 
This will be the subject of a future work and will be discussed elsewhere. 
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Appendix A 

Characteristics of the Function S{x) 

From Eqs. (2-3) and (2-5), the function S{x) is a solution of the following nonlinear 
differential equation: 

A\S(x)+xS'(x)]^ + 9x~-^ = 0. (A-1) 

S[x) 

We can easily see that S{x) = l/3a; is also a solution of Eq. (A-1). This solution is not 

equivalent to S{x) = S{x) since S{0) = (3/4)^/^ from Eq. (2-7). It is veriflable by Eq. (2-7) 

that 

/7r\2/3 
x = x:= i^-j (A-2) 

is a root of the equation S{x) = l/3a;. In other words, the solutions S{x) = l/3a; and 
S{x) = S{x) agree with each other at a; = a; (See Fig. 12). 
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Fig. 12. S{x) and l/3x. 



Let us consider the behaviour of a solution S{x) that is regular in the neighborhood of 
X = 0. Expanding S{x) around x = 0, we have 

S{x) = ao + aix + O(x^). 



Substituting this expression into Eq. (A-1), we have, from the zeroth order of x, 



ao 



1/3 



(A-3) 



(A-4) 



and, from the first order of x, 

ai = -^6-^/1 (A-S) 

Therefore, a solution S{x) that is regular at x = should satisfy 5'(0) = (3/4)^/^. This 
result also means that S{x) is a unique solution that is regular at x = 0. 

Next, we consider a solution S{x) that satisfies S{x) = l/3x. Expanding S{x) around 
X = X, we have 

S = bo + bi{x -x) + h2{x - xf + 0((x - xf). (A-6) 

Substituting this expression into Eq. (A-1), we find, from the zeroth order of (x — x), that 



3x^ 



Then, we have 



S'ix) 



hi 



3x2 



1 
3x 



(A-7) 



(A-S) 



Therefore, the derivative of S{x) at x = x is unique. This result implies that the solutions 
S{x) = l/3x and S{x) = S{x) have the same gradient at x = x. The first order of (x — x) 
in Eq. (A-1) is automatically satisfied. From the second order, we have 



3x^ 



1 27\ 

62 ^ + — =0. 

3x3 IQ 



(A-9) 
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There are two roots 62 = &2± of the above equation, where 

1 1 27 

&2+ = ^^ and 62- = ^^-777- (A- 10) 



This fact means that there are at least two solutions of Eq. (A-1) with the same value and 
the same derivative at a: = a:. The root 62 = ^2+ corresponds to the solution S{x) = l/3a; 
whereas the other root 62 = ^2- corresponds to the solution S{x) = S{x). Since we have 

S + xS' = 2(61 + xb2-){x -x) + 0{{x - xf), (A-U) 

and since 5i < and 62- < 0, we have 

S + xS' >0 for X <x, (A-12) 

S + xS' <0 for x> X. (A-13) 

Eventually, we obtain the differential equation for S{x) as 



S + xS' = —J--3x for x<x, (A-14) 



S + xS' = \/"c~'^^ ^'^^ X > X. (A-15) 

Appendix B 

Dyer-Roeder Equation 



Let us consider a bundle of light rays whose sectional area is given by A. The expansion 
rate 6 oi A along the ray is defined by 

= A, (B-1) 

2A' ^ ' 

where the dot denotes differentiation with respect to an afiine parameter A of the ray bundle. 
Then the evolution of 9 is given by^^-* 

e = -e'-a^- ^R^uP^p", (B-2) 

where p^ is a tangent vector of the light ray, and o"^ is the shear factor defined by 

<y^ = ^v^p.W - li'^.P")'- (B-3) 

Substituting Eq. (B-1) into Eq. (B-2), we have 

^=-(a^ + ^R^uP'-p") VA. (B-4) 
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Let us assume that the hght beam remains far away from clumps and that the contribu- 
tion of a in Eq. (B-4) is neghgible.'^^^"^*^) In addition, we assume that the fraction of matter 
density that is smoothly distributed is given by a, namely, the energy density on the ray is 
given by ap. As mentioned in §1, a is the so-called smoothness parameter. Then, Eq. (B-4) 
reduces to 

VI = -A^Pap^, (B-5) 

where we have used the Einstein equation with the energy momentum tensor (2-2). 
Since the relation between the angular diameter distance and A is given by 

D oc VI, (B-6) 

we obtain 

ij = -inPapD. (B-7) 

When D is given as a function of redshift z, we can rewrite the above equation in the form 
ofEq. (2-19). 

Appendix C 

Derivation of Basic Equations 

From Eq. (2-16), we have 

r='-^. (C4) 



This is identical to Eq. (2-22). From Eq. (2-26), we obtain 

. c 

z ^ 

dD/dz 

as Eq. (2-23). Using C, Eq. (2-19) can be written as Eq. (2-24): 



(C-2) 



(l + zV 
C = -An^-^-^apD. (C-3) 

-"o 
The remaining task is the derivation of Eqs. (2-20) and (2-21) from Eqs. (2-14) and 
(2-18). Differentiating Eq. (2-5) with respect to r, we have 

d,R = -rm-^/\t - t^f^ {S - 2xS') - 
3 f 

3 r 

+ m^'\t-t^f/''S. (C-4) 
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Substituting Eqs. (C-4) and (2-16) into Eq. (2-14)Cwe have 

aril + hi-B + cA; + (i = 0, 



where 



a = -rm'^/^{t - te)'/^ {S - 2xS') 



-rm 



'^Ht-t^r'/'is+xs') 



c = rm ^'\t-tB)*^^S', 



d 



-1/3, 
1+Z 



Hn 



Vl - kr^ 



m 



1/3 



(t - t^r/'s 



Using Eqs. (2-13) and (2-16), we can rewrite Eq. (2-18) as 

2 



-dr-R- 



c 



+ dtdrR 



l + z 



0. 



HodD/dz ' """'" V Ho 
Differentiating Eq. (C-4) with respect to t, we have 

o o'' r 

D o'^ r 

+ -rm-'/\t - t^y/' (25' + xS") - 
3 r 

+ lm'^'{t-tB)-'/'{S + xS% 
where we have used the following equation for S{x) (see Ref.41)): 



x{2SS" + S'^) + 5SS' 
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0. 



Substituting Eqs. (2-14) and (C-11) into Eq. (C-10), we have 

em + /te + gk + h = 0, 



where 



e = rm 



-2/3 



(t - te) 



-1/3 



l + z 



h 



[l + z)dD/dz 



H^ 



[Am'/\t-tB)-''''S\S + xS')\ 



(C-5) 

(C-6) 

(C-7) 
(C-8) 
(C-9) 



(C-10) 



(C-11) 

(C-12) 

(C-13) 

(C-14) 
(C-15) 
(C-16) 

(C-17) 
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If tB{r) is given, from Eqs. (C-5) and (C-13), we have 



where 



m = Fm{m,k,r,z,C) ■= , (C-18) 

ag — ce 

k = F,{m,k,r,z,C)--=^^^^, (C-19) 

ag — ce 



-no 

h = h + l^^±^. (C.21) 

-no 

Appendix D 

{d^R)b, {d^drR% and {dtdlR\ 



We show how to calculate (9^-R, d^drR and dfd^R on the null geodesic that corresponds to 
the observed light ray. Let us consider the differentiation of the null condition (2-14) along 
the null geodesic (t,r) = (tb,Tb) with respect to A as follows: 



d_ 
dX 

Using t = —r, we have 



(drR - \J\ - krA = 0. (D-1) 



{d'^R% = {mR)^-i-^4=^] > P-2) 



2Vl - A;r^ 



b 



where dtdrR is given by Eq. (C-11). 

The expression for d^drR is given by differentiating Eq. (2-3) with respect to t and r as 

where drR is given by Eq. (C-4). 

Finally, {dtdlR\ is given by differentiating the null geodesic equation (2-18) with respect 
to the affine parameter A. Namely, 

— {drRi + dtdrRi^) = 0. (D-4) 

aX 

The above equation leads to 

{dtd^rR)b = {dtdrR)h + tb ^ {drRi' + SdtdrRti - dlRti) ^ , (D-5) 
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where 



t = -^, (D-6) 

' ^ (D-7) 



HodD/dz' 
... _ Atx{1 + zfapD Cd^D/dz^ 
^~ HldD/dz ^ Ho{dD/dzf' ^ '' 

Appendix E 

Initial Conditions at the Center without Distance Modification 



From the null condition, we have dr-R ~ 1 near the center. Thus we must have i? ~ r in 
order that the metric (2-1) is regular at the center. Therefore, from Eq. (2-4), M should be 
given by 

M{r) = '^r'' + 0{r^) (E-1) 

for the finiteness of p at the center. Thus, we have 

m(r) =mo + 0(r). (E-2) 

Hereafter, the subscript describes the initial value at the center. Obviously, ro = -Zo = 0. 
The remaining initial values are mo, A^o and C^q. In addition, the initial value of to should also 
be determined. 
Since 

R^r (E-3) 

near the center, Eqs. (2-3), (2-4) and (2-5) are given by 

(E-4) 

(E5) 
(E6) 



(E-7) 





«7rpo = mo, 


\R ) 


- Hq - /Co + g 

r=0,t=to 


'"'^^'''^ - 8npo' 


where 


Equation (E-5) can be rewritten as 





1 = f^ko + ^Mo, (E-8) 

25 



where 

n^o = -ko/Hl (E-9) 

^Mo = 8npo/{3H^). (E-10) 

Once we fix the value of i7mo at the symmetry center, i7ko, or equivalently ko, is fixed by 
Eq. (E-8). rriQ and po are related to O^q by Eqs. (E-4) and (E-10). tg is given as the solution 
of Eq. (E-6). The remaining initial value is Co = zodD/dz\z=o = zq/ Hq. Differentiating Eq. 
(2-3) with respect to r, we have 

The above equation becomes 

dtdrR = Ho (E-12) 

at r = and t = to- Using this equation, from Eqs. (2-16) and (2-18), we obtain 



1 dfdrR 



Z0 = TT 



Hq OrR 

Thus, we have 



= 1. (E-13) 

r=0,t=to 



Co — Zq — — 

az 



_ zo_ _ J_ 

=0 -^0 Hq 



(E-14) 
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